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In the asymptotic limit Q 2 3> m 2 , the heavy flavor Wilson coefficients for deep-inelastic scattering factorize into 
the massless Wilson coefficients and the universal heavy flavor operator matrix elements resulting from light-cone 
expansion. In this way, one can calculate all but the power corrections in (rn 2 /Q 2 ) k ,k > 0. The heavy flavor 
operator matrix elements are known to NLO. We present the last 2-loop result missing in the unpolarized case 
for the renormalization at 3-loops and first 3-loop results for terms proportional to the color factor Tp in Mellin- 
space. In this calculation, the corresponding parts of the NNLO anomalous dimensions [1,2] are obtained as well. 



1. Introduction 

The unpolarized deep-inelastic double differen- 
tial scattering cross-section can be expressed in 
terms of the structure functions F^cc, Q 2 ) and 
Fl(x,Q 2 ) in case of single photon exchange. In 
the small x region, both structure functions con- 
tain large cc-contributions of the order of 20- 
40 %. Precision extractions of parton distribution 
functions and the measurement of Aqcd there- 
fore require to extend the description of these 
contributions to the O(a^) terms as reached in 
the massless case. The complete NLO corrections 
were calculated semi-analytically in x-space, [3], 
for which a fast implementation in Mellin N- 
space was given in [4]. For Q 2 > 10 m 2 , one ob- 
serves that F2 C (x,Q 2 ) is very well described by 
its asymptotic expression in the limit Q 2 ^> to 2 , 
[5]. In this kinematic range, one can calculate 
the heavy flavor Wilson coefficients, the pertur- 
bative part of the structure functions F2 C (x,Q 2 ) 
and F£ c (x, Q 2 ), analytically, which has been done 
for F| c (a;,Q 2 ) to 2-loop order in [5,6] and for 
Fl c (x,Q 2 ) to 3-loop order in [7]. First steps 
towards an asymptotic 3-loop calculation for 
F| c (x, Q 2 ) have been made by the present au- 
thors by calculating the O(e) terms of the 2-loop 
heavy operator matrix elements (OMEs), [8,9], 
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contributing to the 3-loop heavy flavor Wilson 
coefficients via renormalization. In the present 
paper we also report on new results contributing 
to the 3-loop operator matrix elements and the 
progress towards a full 3-loop calculation of mo- 
ments of the heavy flavor Wilson coefficients. As 
part of the calculation, the Tj.-terms of the even 
moments N — 2. ..10 of the NNLO non-singlet 
(NS) and pure-singlet (PS) anomalous dimensions 
given in [1,2], cf. also [10] in the NS-case, are con- 
firmed in an independent calculation. 

2. Heavy flavor Wilson coefficients in the 
asymptotic limit 

As outlined in Ref. [5], the heavy flavor Wilson 
coefficients 2 in the limit Q 2 ^S> m 2 are obtained 
as a convolution of the light flavor Wilson co- 
efficients with the corresponding massive oper- 
ator matrix elements of the flavor decomposed 
quarkonic and gluonic operators between mass- 
less parton states. Here we consider the level of 
twist-2 operators. The light Wilson coefficients 
are known up to three loops [11] and carry all the 
process dependence, whereas the OMEs are uni- 
versal, process-independent objects. Their log- 
arithmic contributions in TO 2 //i 2 , as well as all 
pole terms in 1/e, are completely determined by 
renormalization. The latter provide checks on the 
calculation. Here, the single pole terms contain 
the respective contributions of the 3-loop anoma- 

2 We consider extrinsic charm production only. 



2 



lous dimensions. 

3. Renormalization 

For the calculation of the bare heavy flavor OMEs 
we use dimensional regularization in D = 4 + e 
dimensions and work in the MS-scheme, if not 
stated otherwise, using Feynman-gauge. The 
renormalization is performed in four steps. For 
mass renormalization we use the on-shell scheme 
[12], whereas charge renormalization is done us- 
ing the MS scheme. For the latter we make the 
requirement that the heavy quark loop contri- 
butions to the gluon self-energy, n(p 2 ,m 2 ), are 
renormalized in such a way that 11(0, to 2 ) = 0, 
cf. [5,6,8]. The origin of the remaining diver- 
gences is twofold. The UV singularities are renor- 
malized via the operator Z-factors, whereas the 
collincar singularities are removed via mass fac- 
torization through the transition functions T. We 
denote the completely unrenormalized OMEs by 

a double-hat, A, and those for which mass and 
coupling renormalization have already been per- 
formed by a single hat. Operator renormalization 
and mass factorization then proceeds via 



A = Z AT 



(1) 



Note that in the singlet case, Eq. (1) should be 
read as a matrix equation, contrary to the NS- 
case. The Z-factors are related to the anomalous 
dimensions of the twist-2 operators via 



7 = /j,dlnZ((i)/dn , 



(2) 



which allows to express the Z-factors in terms of 
the anomalous dimensions up to an arbitrary or- 
der in the strong coupling constant a s := a s / (4ir) 
(cf. [8] up to 0(aj?)). If all quarks were mass- 
less, the transition functions T would be given by 
r = Z _1 , but since we are dealing with diagrams 
containing at least one heavy quark, this equa- 
tion has to be modified in such a way that mass 
factorization is applied to those parts of the dia- 
grams containing massless lines only. Finally, the 
PS and NS terms we are calculating are related 

by 



NS 



Z qq + Z qq 



Jqq 



(3) 



From Eqs. (1,2) one thus can infer that for op- 
erator renormalization and mass factorization at 
0(a,g), the anomalous dimensions up to NNLO, 
[1,2], together with the 1-loop heavy OMEs up 
to 0(e 2 ) and the 2-loop heavy OMEs up to 0(e) 
are needed. The last two quantities enter since 
they multiply Z— and T-factors containing poles 
in e. This has been worked out in some detail in 
Rcf. [8], where we presented the 0(e) terms a^, 

a^g S and Oq^ PS in the unpolarized case. The 

term a gg Q was given in [9]. The missing 2-loop 

0(e) term corresponds to the heavy OME A k 'q 
and was calculated for the first time in Rcf. [13]. 
It contributes through operator mixing to the Tp- 

(3) PS 

term of AW' , which we consider in this paper. 



4. A 



(2) 

gi,Q 



The term A gq ^Q emerges for the first time at 
0(a 2 s ). By applying Eq. (1), one obtains at 0(a 2 s ) 
the renormalized OME 

A (2) i(2) , 7-1,(2) , f 7-1,(1) , \ r -i,(i) 
A gq,Q- A gq + ^gq + {fgg + A gg,Q ) 1 gq 



Here, the term A^Jq, cf. [13], enters through 
mixing. Note that since we consider only terms 
involving at least one heavy quark, we adopt the 
following definition for the anomalous dimensions 



7 = 7(71/ + 1) - i(n f ) 



(4) 



in order to obtain the correct color projection. 
Now we can predict the structure of the unrenor- 
malized result to be 



i (2) 
A gq,Q 



/to 2 \ f 

VTfl) 



_9 Iqq ' 



igq 
2e 



,( 2 ) 



+sa 



(2) 
gq,Q 



(5) 



where the LO and N LO anomalous dimensions are 
given by 

N 2 + N + 2 



lgq F (N -l)N(N + 1) ' 

7 « = C fT J-™ N2+N + 2 Sl 
lgq \ 3 (N — 1)N(N + 1) 



(6) 
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+- 



32 87V 3 + 137V 2 + 27N + 16 
(N - l)N(N + l) 2 



(7) 



and /?o,q = — 4/3Tp . The calculation in 
Mcllin-space in terms of Feynman-parameters is 
straightforward, cf. [6], and a representation in 
Euler-r functions can be obtained even to all or- 
ders in e We reproduce the pole terms of Eq. (5) 
and obtain for the constant term in e 



(2) 

l gg,Q 



T F C F { - 



4 N 2 + N + 2 



S 2 + Sf 



+- 



3 (N - 1)N(N + 1) 

8 8A 3 + 13iV 2 + 277V + 16 

9 (iV- 1)N(N+ l) 2 



Si 
(8) 



27 (TV - l)iV(iV + l) 3 J ' 

Pi = 43A 4 + 105A 3 + 2247V 2 + 230A + 86 

The last term has first been calculated in [13], 
with which we agree. The 0{e) term is new and 
we obtain 



a 



(2) rp r I 2 jV 2 +iV + 2 

F F \ 9(N-1)N(N + 1) 



(-2S 3 



-35 2 S 1 -5 1 3 + 4C 3 -6C25 1 
2 8N 3 + 13N 2 + 27A + 16 



+ 9 (iV - 1)N(N + l) 2 



TT^l 2C2 



+5 2 + 5 2 



4PiSi 



27(7V- 1)A(A + 1) 3 



+ 



4A 



(9) 



81 (iV- 1)A(A + 1) 4 j ' 

P 2 = 248A^ 5 + 8637V 4 + 1927iV 3 + 2582N 2 
+ 1820^ + 496 . 

Eqs. (7-9) are given in terms of harmonic sums, 
[14,15], the argument of which we have set equal 
to N. The representation in Mellin-space allowed 
us to use various analytic and algebraic relations 
between harmonic sums, [16,17], to obtain a more 
compact result. Together with the result of Eq. 
(9), all 2-loop 0(e) terms of the heavy OMEs in 
the unpolarized case are known by now. As a 
last remark, note that we consider charm quark 
contributions here, while for heavier quarks de- 
coupling [18] has to be applied. 



5. Fixed values of N at three loops 

The diagrams we need to calculate are of the 3- 
loop self energy type with the external particle 
being massless and on-shell and contain two in- 
ner scales. One is set by the mass of the heavy 
quark, the other by the Mellin-variable N of the 
operator insertions emerging in the light-cone ex- 
pansion. At this order, new operator vertices ap- 
pear with three and four gluonic lines, for which 
the Feynman-rules had not yet been derived. The 
necessary diagrams are generated using QGRAF 
[19]. The number of diagrams contributing to 

(31 

Ag g , e.g., is 1478 diagrams with one and 489 di- 
agrams with two quark loops. 
The steps for the calculation of these diagrams 
are the following: The diagrams are genuinely 
given as tensor integrals due to the operators con- 
tracted with the light-cone vector A, A 2 = 0. 
The idea is, to first undo this contraction and to 
develop a projector, which, applied to the tensor 
integrals, provides the results for the diagrams for 
the specific (even) Mellin N under consideration. 
So far, we implemented the projector for the first 
5 contributing Mellin moments, N = 2, ...,10, 
where the color factors are calculated using [20]. 
A generalization to higher moments is straightfor- 
ward, however, one quickly runs into computing 
time problems. The diagrams are then translated 
into a form, which is suitable for the program 
MAT AD [21], which does the expansion in e for the 
corresponding massive three-loop tadpole-type 
diagrams. We have implemented these steps into 
a FORM-program, cf. [22], and tested it against 
various two-loop results, including the result for 
A gq Q, Eq. (5), and found agreement. 



6. First results 

The first 3-loop objects being investigated are the 
terms oc T| of the OMEs A^ Q and A^ q . Note 
that in the NS-case, the ± and v terms do not dif- 
fer. All diagrams contain two inner quark loops, 
where the quark to which the operator insertion 
couples is heavy and the other one may be heavy 
or light. The latter two cases can be distinguished 
by a factor n/, denoting the number of light fla- 
vors, in the result. 
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Table 1 

Mcllin moments 2 to 10 for the Tp terms of the O(e ) terms of the unrenormalized 

s(3),PS 

OME A Qq as obtained from MATAD. 

AT ( 3 ). PS I 
_^ a Qq \l% 

36880 736~ 4096 / 76408 U2~ 896 \ 

. 2879939 1118 15488 ^ / 474827503 851 3388 \ 

1 ~ 5467500 ~ 2025 ~ 2025 + Uf y 109350000 ~ 20250 ^ 2 + 2025 j 

146092097 7592 61952 / 82616977 16778 , 1936 ^ \ 

6 C2 C3 + n f \ C2 H C3 

1093955625 99225 s 19845 s } \ 45378900 694575 s 2835 s J 

48402207241 , 1229 43808 / 16194572439593 343781 1369 

272211166080 + 142884 ^ 2 ~ 25515 ^ + Uf \~ 15122842560000 ~ 14288400 ^ 2 + 3645 ^ 

|() 430570223624411 , 319072 802816 

2780024890190625 + 11026125 ^ 2 ~ 735075 Cs 

/ 454721266324013 547424 175616 \ 

+n/ ^~ 624087220246875 ~ 24257475 ^ + 735075 j 



Table 2 

(2) 

Mcllin moments 2 to 10 of 7p S |t 2 as obtained from MATAD and Eq. (11). 



N 


2 


4 


6 


8 


10 


(2)i 
7p S \T% 


5024 
243 


618673 
151875 


126223052 
72930375 


13131081443 
13502538000 


265847305072 
420260754375 



6 1 A (3),PS 

°- L - Qq 

From Eq. (1), we obtain the pole structure of the 
completely unrenormalized PS OME. Considering 
only the Tp terms, one finds 



s(3),PS 
l Q<z 



/^x3e/2 f nf+4 , .(0) (0) 



(«/ + l)^o,Q7^) 



_ j_/2j-n/ (0) (i) 
£ 2 \^ g iqg igq 

l/n/ + U(2) 

£ 



(10) 



+ ~(0) fl (2) \, (3) 



(3), PS 



T 2 - 



where we have written the n / dependence explic- 

A (2) 

itly and with 7 PS being the term oc n"j of the 



(2) 

NNLO anomalous dimension 7 ps \ It is not pos- 
sible to factor out n/ + 1, not even in the triple 
pole term. This is due to the interplay of the pre- 
scription for coupling constant renormalization 
we have adopted, cf. [5,8], and the fact that the 
transition functions T apply to sub graphs con- 
taining massless lines only. We have calculated 
the above term using MATAD for N = 2, 10 and 
all pole terms agree with Eq. (11). In Table 1, we 
show the constant terms in e we have obtained. 

Using Eqs. (8,11), one can obtain moments for 
the 3-loop anomalous dimension 7 ps | T 2 , which 
we show in Table 2. These latter results agree 
with the results from Refs. [2]. Here one has to 
make the replacement n/ — > rif(2Tp), with 
Tp = 1/2, and multiply with 2, to account for the 
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Table 3 

Mellin moments 2 to 10 for the Tp terms of the O(e ) terms of the unrenormalized OME 

s(3),NS 

A qq Q as obtained from MATAD. 

at (3),NSi 

28736 512 2048 / 100096 256 896 

2 — 



! 2048 / 100096 256 896 , \ 

<2 - — C3 + n, (^—^ ~ ^7-C2 + — Caj 



2187 81 

151928299 26542 20096 / 1006358899 13271 8792 \ 

4 5467500 2025 ^ 2 ~ 405 ^ + Hf V 10935000 2025 ^ + 405 

„ 26884517771 1712476, 181504, / 524427335513 856238, 11344 \ 

6 C2 C3 +n f [ C2 H C3 

729303750 99225 s 2835 s ; V 4375822500 99225 s 405 s / 

740566685766263 36241943 632512 / 4763338626853463 36241943 39532 \ 

~ 17013197880000 ~ 1786050 ^ 2 ~ 8505 ^ 3 + Uf V 34026395760000 ~ 3572100 ^ 2 + 1215 
6080478350275977191 2451995507 1543040 

1 124545115080540000 ~ 108056025 ^ 2 ~ 18711 ^ 3 

/ 38817494524177585991 2451995507 96440 \ 

+nf V 249090230161080000 ~ 216112050 + 2673 



Table 4 

Mcllin moments 2 to 10 of 7^ s |t 2 as obtained from MATAD and Eq. (1). 



N 


2 


4 


6 


8 




10 


(2) 
'NS \ T§, 


1792 


384277 


160695142 


38920977797 


27995 


901056887 


243 


30375 


10418625 


2250423000 


1497 


356506500 



different convention for the Z-factors we adopted, 
see Eq. (2). Note that in Tables (1-4), there is 
an overall factor CpTp, which we do not show 
explicitly. 

As an example consider the renormalizcd result 
for the second moment. Applying Eq. (1), we 
obtain 




As in Eq. (11), we observe for all moments in 
the NS and PS case that the terms cx ( 2 disap- 
pear after renormalization, since the correspond- 
ing terms in the light flavor Wilson coefficients 
do not contain even ^-values. This provides us 
with a further check on our calculation, since it 
is a general observation made in many D = 4 
calculations. 

4 (3),NS 

6-2. A qq Q 

For the T|-terms of the heavy OME A ( qq ] ^ S , a 
formula similar to Eq. (11) can be derived, which 
we do not show here. Using again MATAD, we have 
calculated the first 5 non-vanishing moments of 
the completely unrenormalized expression. We 
list the constant terms in e, a q 3 q Q S , of our results 
in Table 3. The pole terms we obtain agree with 
what one expects from Eq. (1) and after renor- 



G 



malization, we again observe that there are no 
C2's left anymore. The values for the moments of 
the terms oc Tp in 7^ 2 g we thus obtain are shown 
in Table 4. These values agree with those in Refs. 
[1,10]. 

7. Conclusions and Outlook 

We calculated the last missing 0(e) contribution 
to the unpolarizcd heavy OMEs for general Mcllin 
variable N at 0(a 2 s ), needed for the rcnormaliza- 
tion at O(etg). Furthermore, we installed a pro- 
gram chain to calculate the corresponding 3 loop 
diagrams to 0(a^) using MAT AD. This chain is now 
existing and we expect first complete results in 
the near future. As a first step, we presented mo- 
ments of the terms oc Tp of the heavy OMEs 

AqqQ S an< ^ ^Qg' PS ' f° r which we found agree- 
ment with the general pole structure expected 
from rcnormalization. This provides us with a 
good check on the method we apply for our cal- 
culation. For the calculation of high moments we 
will apply TFORM, [23], in the future. 
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